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Svetitsky and Yaffe have argued that — if the deconfinement phase transition of a (d+l)-dimensional Yang-Mills 
theory with gauge group G is second order — it should be in the universality class of a d- dimensional scalar model 
symmetric under the center C(G) of G. These arguments have been investigated numerically only considering 
Yang-Mills theory with gauge symmetry in the G = SU(N) branch, where C(G) = Z(iV). The symplectic groups 
Sp(N) provide another extension of SU (2) = Sp(l) to general N and they all have the same center Z(2). Hence, in 
contrast to the SU(N) case, Sp(N) Yang-Mills theory allows to study the relevance of the group size on the order 
of the deconfinement phase transition keeping the available universality class fixed. Using lattice simulations, we 
present numerical results for the deconfinement phase transition in Sp(2) and Sp(3) Yang-Mills theories both in 
(2 + l)d and (3 + l)d. We then make a conjecture on the order of the deconfinement phase transition in Yang- Mills 
theories with general Lie groups SU(N), SO{N), Sp(N) and with exceptional groups G{2), F(4), E(6),E(7), E(8). 
Numerical results for G(2) Yang-Mills theory at finite temperature in (3 + l)d are also presented. 



1. Introduction and Overview 

The Yang-Mills theory with gauge group G has 
a finite-temperature phase transition where the 
confined colorless glueballs deconfine in a gluon 
plasma. This transition is signalled by the spon- 
taneous breaking of a global symmetry related to 
the center G(G) of G. The corresponding order 
parameter is the Polyakov loop which transforms 
non-trivially Q{x)' — z$(x) under a global cen- 
ter transformation characterized by the center el- 
ement z E C(G). The expectation value of the 
Polyakov loop ($) = exp(— (3F) is related to the 
free energy F of a static quark as a function of 
the inverse temperature (3 — 1/T. In the low- 
temperature confined phase the center symmetry 
is unbroken, i.e. ($) = 0, and hence the free en- 
ergy of a single static quark is infinite. In the 
high-temperature deconfined phase, on the other 
hand, the center symmetry is spontaneously bro- 
ken, i.e. ($) 7^ 0, and the free energy of a quark 
is finite. Integrating out the spatial degrees of 
freedom of the (d + l)-dimensional Yang-Mills 
theory, one can write down an effective action 
for $. It describes a scalar model with global 
symmetry G(G) in d dimensions. The G(G)- 
symmetric confined phase of the gauge theory 
corresponds to the disordered phase of the scalar 



model, while the C(G)-broken deconfined phase 
has its counterpart in the ordered phase. Svetit- 
sky and Yaffe P argued that the interactions in 
the effective description are short ranged. Hence, 
if the deconfinement phase transition is second 
order, approaching criticality, the details of the 
underlying short-range dynamics become irrele- 
vant and only the center symmetry C (G) and the 
dimensionality d of space determine the universal- 
ity class. Thus, one can exploit the universality of 
the critical behavior to use a simple scalar model 
to obtain information about the much more com- 
plicated Yang-Mills theory. On the other hand, 
if the phase transition is first order, the corre- 
lation length does not diverge and there are no 
universal features. In this case the G-symmetric 
Yang-Mills theory in (d + 1) dimensions and the 
G(G)-symmetric d-dimensional scalar model do 
not share the same critical behaviour. 

For example, the (d + l)-dimensional SU(2) 
Yang-Mills theory has center symmetry Z(2) and 
hence the effective theory is a (/-dimensional Z(2)- 
symmctric scalar field theory for the real-valued 
Polyakov loop. The simplest theory in this class 
is a $ 4 theory with the Euclidean action 
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The scalar potential is given by V($>) = a$ 2 + 
fe$ 4 , where b > for stability reasons. Indeed, for 
a = this theory has a second order phase transi- 
tion in the universality class of the d-dimensional 
Ising model. However, this does not guaran- 
tee that the deconfinement phase transition in 
SU(2) Yang-Mills theory is second order. In par- 
ticular, one could imagine an effective potential 
V($) = a$ 2 + 6$ 4 + c$ 6 , involving also a $ 6 
term which is marginally relevant in three dimen- 
sions. The coefficient c has to be positive in order 
to ensure that the potential is bounded from be- 
low, but the coefficient b of the $ 4 term can now 
become negative. Then the phase transition be- 
comes first order. The above argument extends 
straightforwardly to a generic gauge symmetry 
group G with center C(G). Hence the order of 
the deconfinement transition is a dynamical issue 
that can be addressed only by lattice simulations. 

The Yang-Mills theory on the lattice is natu- 
rally formulated in terms of group elements while 
in the continuum the fundamental field is the 
gauge potential, living in the algebra. An alge- 
bra can generate different groups, however it is 
natural to expect that lattice Yang-Mills theo- 
ries whose gauge groups correspond to the same 
algebra have the same continuum limit. Hence, 
instead of SO(N) we consider its covering group 
Spin(N). Keeping this in mind, we look at the 
center subgroups C(G) of the various simple Lie 
groups G 

C(SU(N))=Z{N); C(Sp(N)) = Z(2) (2) 

fZ(2); TV odd 

C(SO{N))^C(Spin(N)) =1 Z(2) 2 ; N — Ak (3) 

[z(4); Y=4fc+2 

C(G(2)) = C(F(4)) = C(E(8)) = {1} (4) 
C(25(6)) = Z(3); C(E(7)) = Z(2) (5) 
Numerical simulations in (2 + 1) and (3 + 1) di- 
mensions have been performed for SU(N) Yang- 
Mills theory in order to investigate the order of 
the deconfinement transition and - in case it is 
second order - to check the validity of the conjec- 
ture of Svetitsky and Yaffe. The currently known 
results are: 

• (3 + 1) dimensions. The SU(2) Yang-Mills the- 
ory has a second order deconfinement transition. 
Consistent with the Svetitsky- Yaffe conjecture, it 



is in the universality class of the 3d Ising model. 
SU(N) Yang-Mills theories with N = 3,4,6,8 
have a first order deconfinement phase transition: 
hence there are no universal features. 
• (2 + 1) dimensions. Lowering the dimensional- 
ity of space makes the fluctuations stronger. In- 
deed SU(N) Yang-Mills theories with N = 2, 3 
and, perhaps, also Af = 4 | have a second 
order deconfinement phase transition. Also in 
these cases the Svetitsky- Yaffe conjecture is veri- 
fied and one finds that the 2d universality classes 
are, respectively, those of the Ising, 3-state Potts, 
and Ashkin- Teller models. 

However the branch of SU(N) groups is not 
a good choice to study the relation between the 
order of the deconfinement phase transition and 
the size of the group. In fact, in this case, when 
increasing the size (N 2 — 1) of the group also 
the center Z(iV) changes. In order to disentangle 
these two features we have considered the Yang- 
Mills theory with gauge group Sp(N). Since all 
Sp(N) groups have the same center Z(2), we now 
keep fixed the available universality class and we 
can directly study the relevance of the size of the 
group on the order of the deconfinement transi- 
tion. Finally, this also allows us to explore the 
deconfinement transition in Yang-Mills theories 
with a gauge symmetry different from SU(N). 
The results presented in these proceedings are 
published in |3j 

2. The Symplectic Group Sp(N) 

The group Sp(N) is the subgroup of SU(2N) 
which leaves the skew-symmetric matrix 

J=(\ I ) = i°2 ® a, (6) 

invariant. Here 02 is the imaginary Pauli matrix 
and 1 is the N x N unit-matrix. The elements 
U e SU(2N) belonging to Sp(N) satisfy the con- 
straint 

U* = JUrf. (7) 

Hence U and U* are related by the unitary trans- 
formation J: consequently, the 2iV-dimensional 
fundamental representation of Sp(N) is pseudo- 
real. The matrix J itself also belongs to Sp(N), 
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implying that in Sp(N) Yang-Mills theory charge 
conjugation is just a global gauge transformation. 
This property is familiar from SU(2) = Sp(l) 
Yang-Mills theory. 

The constraint eq. (JJJ implies the following form 
of a generic Sp(N) matrix 



U 



W 

-X* 



X 

w* 



(8) 



where W and X are complex N x N matrices 
such that WW^+XX* = 1 and WX T = XW T . 
Since center elements are multiples of the unit- 
matrix, eq.® implies W = W* oc 1: hence, the 
center of Sp(N) is Z(2). 

Writing U = exp(iif), where H is a Hermitean 
traceless matrix, eq.{7J) implies that the genera- 
tors H of Sp(N) satisfy the constraint 

H* = -JHJ^ = JHJ. (9) 

This relation leads to the following generic form, 



H 



A 

B* 



B 

-A* 



(10) 



where A and B are N x N matrices such that 
A = A^ and B = B T . Since A and B have, 
respectively, N 2 and (N+1)N degrees of freedom, 
the dimension of the group Sp(N) is (2N + l)N. 
There are N generators of Sp(N) that can be 
simultaneously diagonalized and so the rank is 
N. The N = 1 case is equivalent to SU(2), while 
the N = 2 case is equivalent to SO (5), or more 
precisely to its universal covering group Spin(5). 

3. Sp(N) Yang-Mills Theory on the Lattice 

The construction of Sp(N) Yang-Mills theory 
on the lattice is straightforward. The links U x „ € 
Sp(N) are group elements in the fundamental 
{27V} representation. We consider the standard 
Wilson plaquette action 



S[U] 



• y p 



(11) 



where g is the bare gauge coupling and Up = 
Ux.fJJx+^vUl^ ^Ul v is the plaquette. The par- 
tition function then takes the form 



VUexp(-S[U]), 



(12) 



where the path integral measure is the product of 
the Sp(N) Haar measures of each link. Both the 
action and the measure are invariant under gauge 
transformations 

= n x u x ^ni + ~, n x G s P (N) (13) 

The Polyakov loop is defined by 

N t 

<P S = Tr(vY[U SAd+1 ) (14) 
t=i 

where Nt = 1/Tis the extent of the lattice in Eu- 
clidean time, which determines the temperature 
T in lattice units. The lattice action is invariant 
under global Z(2) center symmetry transforma- 
tions U' S Nt d+1 = -Ua t N t ,d+i, while the Polyakov 
loop changes sign, i.e. $^ = — As a conse- 
quence of the center symmetry, the expectation 
value of the Polyakov loop always vanishes, i.e. 
(<fr) = 0, even in the deconfined phase. This is 
simply because spontaneous symmetry breaking 
does not occur in a finite volume. However, the 
Polyakov loop probability distribution does 
indeed allow one to distinguish confined from de- 
confined phases. In the confined phase p($) has a 
single maximum at $ = 0, while in the deconfined 
phase it has two degenerate maxima at $ ^ 0. If 
the deconfinement phase transition is first order, 
the confined and the two deconfined phases co- 
exist and one can simultaneously observe three 
maxima close to the phase transition. At a sec- 
ond order phase transition, on the other hand, the 
high- and low-temperature phases become indis- 
tinguishable. The two maxima of the deconfined 
phase merge and smoothly turn into the single 
maximum of the confined phase. Three coexist- 
ing maxima then do not occur in a large volume. 

As physical quantities useful in the finite-size 
scaling analysis presented below, we also intro- 
duce the Polyakov loop susceptibility 



X 



as well as the Binder cumulant 
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The susceptibility measures the strength of fluc- 
tuations in the order parameter while the Binder 
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cumulant measures the deviation from a Gaus- 
sian distribution of those fluctuations. We also 
consider the specific heat which takes the form 

1 



C v = 



L d Ni 



as 2 ) - (syy 



(17) 



In a finite volume Cy has a maximum close to 
the critical coupling AN/g 2 of the infinite volume 
theory. We denote the value of the specific heat 
at the maximum by Cy ax . Another interesting 
observable is the latent heat 

L H =j^-{{S) c -{S) d ), (18) 

which measures the difference of the expectation 
values of the action in the confined and the de- 
confined phase. In the large volume limit Lh and 
are related by (SJ 



r^raax 



T 2 

C max = L d N H 



4. Numerical Results 
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The Sp(N) lattice Yang-Mills theory with the 
standard Wilson action can be simulated us- 
ing the Cabibbo-Marinari method 6 of alter- 
nating heat-bath [7] and microcanonical overre- 
laxation |8I9I10| algorithms in various SU(2) = 
Sp(l) subgroups. In the next two subsections we 
report on the results of numerical simulations in 
Sp(2) and Sp(3) Yang-Mills theories in (2 + l)d 
and (3+l)d. We denote with L and N t the spatial 
and the temporal lattice extensions, respectively. 

4.1. Sp(2) Yang-Mills Theory 

As a first step we have scanned the expecta- 
tion value of the plaquette in order to check if 
the strong and the weak coupling regimes are sep- 
arated by a bulk transition. Fortunately, we find 
no bulk phase transition that might interfere with 
the study of the deconfinement transition both in 
(2 + l)d and (3 + l)d. Our Monte Carlo data for 
the expectation value of the plaquette are com- 
pared with analytic weak and strong coupling ex- 
pansions in figure ^ 

In (2 4- l)d we observe a second order decon- 
finement transition, signalled by the broadening 
of the probability distribution of $ and, hence, by 
the increase of the Polyakov loop susceptibility % 
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Figure 1. Monte Carlo data from hot and cold 
starts for the plaquette in (2 + l)d and (3 + l)d 
Sp(2) Yang-Mills theory compared to analytic re- 
sults in the weak and strong coupling limits. 

at criticality. A finite size scaling analysis con- 
firms the expectation that the universality class 
is that of the 2d Ising model. Figure [5] shows the 
collapse on a single curve of the data for (|$|) col- 
lected on lattices of different sizes L 2 x 2 and at 
various couplings g 2 . The variable x — {g 2 /g 2 — l) 
is a measure of the distance from the critical cou- 
pling g 2 (N t — 2). The critical exponents v, (3 and 
7 are characteristic of a universality class: in this 
case they have been fixed to those of the 2d Ising 
one. In figure we also plot rescaled data for 
5(7(2) Yang-Mills theory in (2 + l)d which has a 
second order deconfinement transition in the 2d 
Ising universality class. The two sets agree excel- 
lently. In (3 + l)d - contrary to what one might 
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Figure 2. (2 + \)d Sp(2): finite-size scaling plot 
for {\^\)L f) l v . SomeSU{2) data are included too. 

have expected - the probability distribution of $ 
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in the critical region clearly shows the coexistence 
of the symmetric and of the broken phases. This 
indicates that the deconfmement transition is first 
order. Figure [3] shows the susceptibility x as a 
function of the gauge coupling for different spatial 
sizes L, keeping N t — 2 fixed. It clearly turns out 
that, at the critical coupling 8/g% = 6.4643(3), x 
scales with the spatial volume L 3 . This quantita- 
tively confirms the first order nature of the phase 
transition. Figure 0] shows the maximum of the 

0.020 I 
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8/g' 

Figure 3. Scaling of x i> n Sp(2) Yang-Mills theory 
onL 3 x2 lattices . We estimate 8/g* = 6.4643(3). 

specific heat per volume, Cy ax /L 3 , as a function 
of the inverse volume 1/L 3 . The linear behavior 
is characteristic of a first order phase transition. 
A linear extrapolation of Cy ax /L 3 to the infinite 
volume limit (see eq.^HJ) is consistent with a di- 
rect measurement of Lh- This again supports the 
first order nature of the transition. 

We have also checked if the deconfmement tran- 
sition for (3 + l)d Sp(2) Yang- Mills theory re- 
mains first order in the continuum limit. To test 
for scaling, we have measured the dimensionless 
ratio T c /^/a between the deconfmement temper- 
ature T c and the square root of the string tension 
a at T = 0. Our data indicate that we are in 
the scaling regime, showing a scaling behaviour 
proportional to the lattice spacing squared. 

4.2. Sp(3) Yang-Mills Theory 

The results of Sp{2) Yang-Mills theory show 
that in (2 + l)d fluctuations are stronger than in 
(3 + l)d and the deconfmement transition is sec- 
ond order. Expecting that the larger the group 
the weaker the fluctuations, we have also in- 
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Figure 4. Approach to the infinite volume limit of 
Cy ax jL 3 for (3+1) d Sp{2) Yang-Mills theory. 
The linear extrapolation is in agreement with the 
measured latent heat Lh- 

vestigatcd the deconfmement transition in 5*15(3) 
Yang-Mills theory. Consistent with this picture, 
we find that in (2 + I) dimensions Sp(3) Yang- 
Mills theory has a first order deconfmement tran- 
sition. The probability distribution of $ in the 
critical region indeed displays the coexistence of 
the broken and of the symmetric phases. Equiv- 
alently, one can say that, since the number of 
colorless glueball states is almost independent 
of the gauge group, the larger number of Sp(3) 
gluons w.r.t. Sp(2) gluons in the deconfmed 
phase, increases the difference between the rel- 
evant degrees of freedom on the two sides of the 
phase transition. This drives deconfmement to 
take place with a discontinuous first order tran- 
sition. In (3 + l)d - similar to the Sp{2) case - 
Sp(3) Yang-Mills theory deconfines with a first 
order transition. Figure EI shows tunneling events 
between the coexisting symmetric and broken 
phases as a function of Monte Carlo time t MC . Fi- 
nally, we note that Sp(3) Yang-Mills theory has no 
bulk phase transition both in (2+l)d and (3+1 )d. 

5. Conjecture 

Our numerical results show that (2 + l)d Sp(2) 
Yang-Mills theory has a second order deconfme- 
ment transition in the 2d Ising universality class. 
However, (3 + l)d Sp{2), (2 + l)d and (3 + l)d 
Sp(3) Yang-Mills theories deconfine with a non- 
universal first order transition. Hence, despite the 
fact that a universality class is available, Yang- 
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Figure5. (3+l)d Sp(3): tunneling events between 
the symmetric and the broken phases. 

Mills theory can have a non-universal first or- 
der deconfinement transition. A non-trivial cen- 
ter plays no role in determining the order of this 
transition. Instead our Sp(N) and the SU (N) re- 
sults indicate that the order of the deconfinement 
transition is a dynamical issue related to the size 
of the gauge group. We conjecture that the dif- 
ference in the number of the relevant degrees of 
freedom between the confined phase (color sin- 
glet glueballs) and the deconfined phase (gluon 
plasma) determines the order of the deconfine- 
ment transition. Thus, we expect that in (3 + l)d 
only SU(2) Yang-Mills theory has a second or- 
der deconfinement transition; in (2 + l)d, due to 
stronger fluctuations, only SU(N), N = 2, 3 and, 
perhaps, N = 4, and Sp(2) Yang-Mills theories 
should have second order transitions. According 
to this picture, E(6) and E(8) Yang-Mills theories 
should also have a first order transition due to the 
large size of the groups: 78 and 133 generators, 
respectively. For Yang-Mills theories with triv- 
ial center gauge groups G(2), F(4), E(8) [TTj the 
large number of generators also suggests the pres- 
ence of a first order transition even if no symme- 
try can break down. Numerical results for G(2) 
Yang-Mills theory in (3 + l)d indeed show the 
presence of a first order finite temperature phase 
transition. In figure we plot the Monte Carlo 
history of the Polyakov loop at finite temperature 
in the critical region: several tunneling events 
clearly show up. The probability distribution has 
a two-peak structure indicating the presence of a 
first order finite temperature transition. 




Figure 6. (3 + l)d G(2): Polyakov loop Monte 
Carlo history with tunneling events between two 
coexisting phases. The data refer to a simulation 
on a 18 3 x 6 lattice. 
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